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Abstract 

The random sequential adsorption (RSA) model is modified to describe damage and crack 
accumulation. The exclusion for object deposition (for damaged region formation) is not for the 
whole object, as in the standard RSA, but only for the initial point (or higher-dimensional defect) 
from which the damaged region or crack initiates. The one-dimensional variant of the model is 
solved exactly. 
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I. INTRODUCTION 



Random sequential adsorption (RSA) has found many applications,- 1 ^ notably in irre- 
versible surface deposition.-^ In this article, we consider a variant of the RSA model mo- 
tivated by study of damage accumulation, e.g., multicracking. The first adaptation of the 
RSA model in mechanical-engineering literature to analyze fiber fragmentation in a uni- 
directional composite,- was followed by defect modeling in thin coatings.-^ 1 ^ In parallel, 
RSA-type models were also utilized in the statistical mechanics literature to explore ap- 
proaches to damage development in several situations.— ^^2. 

Formation of a crack relieves stress in the nearby material, thus creating a region along the 
crack, possibly exclusive of the crack tip/edge area, in which new cracks cannot be initiated. 
Cracks formed in a thin coating can be approximated as a collection of one- dimensional (ID) 
"rods" (segments), subject to the no-overlap restriction as they form sequentially, similar to 
the traditional RSA. Experiments on multiple cracking in metat^^JL and nonmetat^ 1 ^^ 
coatings have been modeled in this way. Furthermore, experiments^ 1 ^ 1 ^ 1 ^ on load-induced 
damage in uniaxial fiber materials have also been interpreted in terms of RSA-type models 
that involve mathematics^ similar to that encountered in the present article. 

Generally, in materials subject to stress, cracks are initiated from static microscopic 
defects.— 1 ^ In the standard formulation of the RSA model, newly formed (or deposited on 
surfaces) objects cannot overlap objects formed earlier. Here we consider a variant of the 
RSA model that is arguably a better representation for damage accumulation in materials. 
Indeed, since the formed cracks relieve stress, except perhaps at their edges (as long as they 
grow), then new cracks cannot be initiated at or near the existing cracks. Therefore, as will 
be described later, the RSA model has to be modified to allow for "exclusion" of only the 
initial region of the crack. We do allow new cracks to be initiated close to the tips of the 
earlier formed cracks. 

Despite the seeming simplicity of the RSA model, exact results are usually available only 
in one dimension^^ and in other special lattice geometries.- In this article, we consider a ID 
model of deposition of segments of length £ on a line, and we allow exclusion to apply only 
to the head points of the segments. We argue that this model describes the accumulation 
of cracks which can reach lengths up to I in ID geometries. 

In Section [Til we describe the model and then in Section HTT| carry out an exact calculation 
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of the open-interval probability, which allows us to evaluate the density of cracks as a 
function of time, t. In Section IIVt we analyze the probability distribution of gaps, exactly 
and asymptotically for large times. 

II. RSA OF SEGMENTS WITH SINGLE-END EXCLUSION 

Our model consists of the deposition of segments of length £ on a linear substrate, as 
illustrated in Figure 0Ja). The segments are transported to the substrate at the rate R 
(number of deposition attempts per unit length per unit time). Each segment has a "head" 
point, and the segment orientation is otherwise random: right-to- left or left-to- right segments 
arrive at the rates R/2 per unit length. Unlike the ordinary RSA, exclusion applies only to 
the adsorption of head points, which can be deposited only in areas that are not covered 
by the previously deposited segments. The lower pairs of segments in Figure [fla) show two 
such possible deposition sequences. However, the segments with head points arriving above 
the covered areas are rejected (their deposition attempts fail and they are transported away 
from the substrate). One such example is illustrated in Figure [T]^a). 

Figure [D(b) illustrates the equivalence of the deposition model to crack formation. The 
head points of the segments correspond to the initiation points of the cracks in a ID structure, 
e.g., fiber .-^^ For tractability (otherwise our model would not be exactly solvable by the 
utilized technique), we assume that the cracks grow up to a fixed length £, except when 
they run into the cracks formed earlier, as shown in the figure for two configurations, each 
corresponding to the first two segments deposited in the corresponding locations in Figure 
[TJ^a). Furthermore, we consider the situation where the cracks grow instantaneously once 
initiated, as compared to the time scales of new crack initiation. 

Since there is a zone near a crack with reduced stress (except perhaps near the tip/edge), 
new cracks cannot initiate nearby. In fact, in real materials higher-dimensional damage 
configurations should be considered. They can originate not only from point-like defects and 
not all be identical in size and shape.— ^ A simple illustration is given in Figure [T](c). Here 
circular damage regions originate from random centers in a plane. Additional damage can 
be initiated only from a center-point located in the previously undamaged areas. Figure [0(c) 
illustrates that the present models can be quite close to the ordinary RSA models. Indeed, 
unless the exclusion point is taken off-center or the objects are made less symmetrical, the 
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model in the figure is nearly identical to the ordinary RSA of half-diameter circles (with no- 
overlap exclusion): the only difference is in the way the "covered" area is counted, whereas 
the center deposition kinetics is identical for both models. 

We note that in the crack-formation nomenclature the parameter R corresponds to the 
rate (per unit time) of crack initiation per unit length of the ID structure. Further discussion 
(and literature citation) is available^ on how such statistical-mechanical parameters relate 
to the actual system properties. In this work we consider the simplest model, the main 
virtue of which is the possibility of deriving analytically exact or asymptotic results. More 
complicated models, which can be realistically compared with experimental data, would 
require large-scale numerical simulations and will be the subject of future work. 

The quantities of interest for the model just defined are as follows. It is interesting to cal- 
culate the fraction of the total length that is undamaged, F(t), which is expected to decrease 
monotonically from F(0) = 1 to F(oo) = 0. A second quantity of interest is the density 
of deposited segment heads (the density of cracks initiated), N(t). Obviously, iV(0) = 
because, for simplicity, we assume that initially the system is undamaged. However, we will 
establish that this quantity actually diverges logarithmically as t — > oo. The calculation of 
F(t) and N(t) will be carried out exactly in the next section. 

Another quantity of interest is the density (per unit length), Q(x,t)dx, of gaps (between 
the damaged areas) of sizes from x to x + dx. The quantity, Q(x,t), gives us the size- 
distribution of undamaged regions in the ID material sample. The total density (number 
per unit length) of gaps between the damaged regions is then 

POO 

G(t) = / Q(x,t)dx. (1) 
Jo 

Note that Q(x, 0) = and, obviously, 

G(t <C (Biy 1 ) ^Rt. (2) 

In addition to exact results, it is of interest to explore the small-gap statistics leading 
to the derivation of the asymptotic results for Q(x,t) in the large-time regime, when the 
total number of gaps has reached the saturation value G(oo). This matter will be further 
elucidated in Section HVl where the results for the gap density are presented. 
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III. EMPTY-INTERVAL PROBABILITY 



Exact results for the present model can be derived for the empty-interval probability, 
P(x,t). This quantity represents the probability that a randomly chosen interval of length 
x is fully undamaged, i.e. open for deposition of segment-heads. Since such an interval could 
be a part of a longer empty region, P(x,t) is not proportional to the gap density Q(x,t)dx 
introduced earlier. These quantities will be related in Section IIVI However, one can easily 
convince oneself that the fraction of the total length that is undamaged is given by 

F(t) = P(0,t). (3) 

The probability P(x,t) satisfies the equation 

x+i 

dP{x,t) 



Of 



RxP(x,t)-R J dyP(y,t), (4) 



where the first term on the right-hand side accounts for depositions where a head-point 
lands in the interval x. The second term on the right-hand side corresponds to deposition 
of segments with the head-point landing outside the interval x, which requires that a larger 
interval, of length y (up to x + £), be open for deposition. 

For the initial condition P(x, 0) = 1, the x dependence in fl4]) can be eliminated by the 
Anzatz 

P(x, t) = exp{-Rxt) F{t) . (5) 

Substituting in (j3J, one obtains a differential equation for F(t) with the initial condition 
F(0) = 1. The result is 

t 

P(x, t) = exp -Rxt + J ex P(~^) ~ l du _ (g) 

o 

For the fraction of the undamaged area, we have 

F{t) = ^e~*W , (7) 

oo 

where Ei(t) = J - — du is a standard exponential integral function, and 7 = 
0.57721566490 ... is the Euler gamma constant.— For large times, the fraction of the un- 
damaged (available for deposition) length vanishes according to 

e-> I e- m < 



n<»wr) = Wt {i-^ E+ -)~ m , <8) 



where from now on the sign « will denote the large-time asymptotic behavior. 

While the total undamaged length shrinks to zero at large times, the total density of the 
deposited segments (the total density of cracks initiated), N(t), actually diverges logarith- 
mically, 

N (t > {My 1 ) « — \n(R£t) . (9) 

This behavior, characteristic of the present model, is opposed to ordinary ID RSA. Indeed, 

t t 
N(t) = R F(v)dv = R / exp 



Lo 



This can be also written 

m) = - 

where the first term dominates the divergence 



mt 

E ^\n{R£t) - I du—e- u e- E ^ 
J u 





(10) 



(11) 



IV. THE GAP DENSITY DISTRIBUTION 

In order to relate the gap density distribution Q(x,t) to the empty interval probability 
P(x, t), let us discretize the distances in the problem in steps Ax, where ultimately, Ax — > 0. 
Specifically, the linear substrate will become a lattice of spacing Ax, and all the other lengths, 
including x and £, will be assumed multiples of Ax. 

Let us denote by fi(x, t) the probability that a randomly chosen interval of length x = 
nAx has its rightmost lattice site blocked, while all the other, n — 1, lattice sites are empty. 
We can write a relation between probabilities, 



P(x - Ax, t) = P(x, t) + Q(x, t) 



which, in the limit of small Ax, gives 



dP(x,t) 
dx 



Ax. 



(12) 



(13) 



Let us now define the probability $(x, t) that a randomly chosen interval of length x = 
nAx has both its rightmost lattice site and its leftmost lattice site blocked, while all the 
other, n — 2, lattice sites are empty. The probability relation is then 



Vt(x - Ax, t) = Q(x, t) + $(x, t) , 



(14) 



from which we get, for small Ax, 



$(x,t) = ^M(A.) 2 . (15) 



Since the density (per unit length) of fixed-length intervals on the lattice is 1/Ax, for the 
density of gaps of size x we get 



Finally, we get 



Q(x, t)Ax = $(x, t)/Ax = d2P g ^^ Ax . (16) 



Q( x , t ) = = {Rtli)e- mt - E ^ mt ^ . (17) 

ox 1 



The total density of gaps per unit length, cf. ([T]), is thus 

G{t) = e~ E ^ mt) -y I . (18) 

This function is shown in Figure [2j 

Since exact solutions are not possible in most geometries, considerations that can yield 
asymptotic large-time behavior play an important role in the study of RSA models. We 
point out that in the opposite, short-time limit, when deposition events are initially not 
correlated, a low-density approximation leads to relations of the type ([2]). The large-time 
behavior, however, is less straightforward. It can be analyzed to various degrees, depending 
on the model details, by considering the dynamics of gaps small enough such that their 
evolution can be described approximately in term of their size and shape. Specifcially, for 
ordinary RSA one considers^^i the regime of late times when most of the remaining gaps 
are small enough to accommodate at most one depositing object. One further assumes 
that the distribution of these gaps at some large to, from which time on the small-gap- 
dominated behavior sets in, is smooth with respect to size, and their shape distribution is 
also known (since we know the geometry of the deposited objects surrounding the gaps) and 
is similarly smooth. The kinetics of these gaps, blocked by the arriving depositing object, is 
then evaluated2&£L in order to obtain the evolution of the gap size and shape distributions 
for t > to, from which other deposit properties can be found. 

In our case, it is natural to consider the large-time behavior of the gap distribution under 
the assumption that it is dominated by gaps not larger than the segment size I. Indeed, such 
gaps cannot be fragmented by segment deposition: they can only be shortened. Therefore, 
the asymptotic regime is for times t > t , see Figure [21 late enough for the function G(t) to 



have reached its constant value. However, other quantities, including the gap distribution 
Q(x,t), with x < £, still vary with time. The equation that describes this variation is 

dQ(x,t) 



dt 



-RxQ(x,t) + R J Q(y,t)dy. (19) 



Here, not only is the equation already approximate, but we seek the asymptotic large-time 
solution in the regime dominated by small gaps, for which we can further replace I — > oo in 
the upper limit of the integration, as shown, because I actually plays no role in the evolution 
of small gaps (most of the length of the deposited segments covers already blocked area and 
is not active in further blocking). 

The solution of the resulting equation is 

Q(x,t) « qte~ Rxt , (20) 

where q is a constant that cannot be determined from the large-time analysis. It can be 
found by assuming properties for the distribution at a reference time to, as mentioned for 
the standard RSA problem. Alternatively, we can calculate (or numerically estimate) some 
less detailed property related to the gap density distribution, for instance the fraction of the 
total length that remains open for deposition 

oo 

F(t) = j x Q( x ,t)dx*^L- t . (21) 
o 

Here we just use the exact large-time expression ([8]) to get 

q = R e -y£ J (22) 
which, finally, gives the asymptotic form of the gap density distribution, 

Q{x,t) ?z Rte- Rxt ~"</£, (23) 

consistent with the exact result ( flTfl . 

In summary, we have studied an exactly solvable RSA-type model motivated by the 
problem of accumulation of cracks in a ID system. The exclusion condition applies only 
to the head-points of the segments, deposition of which models crack formation. Study of 
more realistic systems in higher dimensions will require numerical simulations. We point 
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out that even in the simplest ID model considered here, some interesting quantities are 
not amenable to exact solution. For example, we can define the probability, B(x,t)dx, 
of the blocked intervals of length from x(> £) to x + dx, similarly to the distribution of 
the unblocked gaps, Q(x,t)dx. The quantity B(x,t) is of particular interest because it 
is expected to have a delta-function singularity at x — £, corresponding to a final (initially 
increasing, but eventually decaying as t —> oo) density of isolated deposited segments forming 
blocked-intervals of length exactly I. We were not able to calculate B(x,t) exactly. Finally, 
we explored the large-time small-gap-dominated asymptotic behavior along the lines of the 
derivation appropriate to standard RSA models. 
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FIG. 1: Definition of the model: (a) Examples of possible configurations of sequential deposition 
of segments with head-point exclusion on a line. The curly arrow marks an arriving segment 
the deposition attempt of which will be rejected, (b) Formation of cracks that maps onto the 
deposition problem corresponding to the deposited pairs of segments in panel (a) above. The 
white area denotes the parts of the later- formed cracks that overlap already formed cracks, (c) A 
higher-dimensional example: Here circular damage areas are assumed to spread from center points. 
New damage area cannot be initiated from a center-point that overlaps earlier damaged areas. A 
configuration that will not be realized is marked by a curly arrow. As further explained in the 
text, this model illustrates that the point-exclusion RSA can be nearly identical to the ordinary 
RSA (here, of deposition of half-radius circles) especially when the objects and the exclusion point 
location are symmetrical. 
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